cd z:\R\concentration\

use weidmann_jcr-2009-53-4.dta, clear

label variable CONC "Concentrated"
label variable SPATCONC "Territorial concentration"
label variable POPCLDIST "Population dispersion"
label variable lr3 "Power balance"
label variable ldistpop "Capital distance"
label variable mtgrp "Mountains"
label variable lgdp96l "GDP p.c."
label variable lpopl1 "Country population"
label variable democracy "Democracy"
label variable peaceyrs "Peace years"
label variable year "Year"

estimates clear

logit ethons05d CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl*, nolog cl(cow)
est store m1

logit ethons05d SPATCONC CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl*, nolog cl(cow)
est store m2

logit ethons05d POPCLDIST CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl*, nolog cl(cow)
est store m3

logit ethons05d SPATCONC POPCLDIST CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl*, nolog cl(cow)
est store m4

* VIF computation using the collin package from http://www.ats.ucla.edu/stat/stata/ado/analysis
collin ethons05d SPATCONC POPCLDIST CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl*

logit ethons05d SPATCONC POPCLDIST lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl* if CONC!=1, nolog cl(cow)
est store m5

outreg2 [*] using table, replace tex alpha(0.01, 0.05) label keep(CONC SPATCONC POPCLDIST lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs)


* Robustness checks

* Only groups with dyadic power share of 0.01 or above
logit ethons05d SPATCONC POPCLDIST CONC lr3 ldistpop mtgrp lgdp96l lpopl1 democracy  year peaceyrs spl* if r3>=0.01, nolog cl(cow)
logit ethons05d SPATCONC POPCLDIST lr3 ldistpop mtgrp lgdp96l lpopl1 democracy year peaceyrs spl* if CONC!=1 & r3>=0.01, nolog cl(cow)

* cross-sectional design
logit s_ethons05 SPATCONC POPCLDIST CONC s_lr3 ldistpop mtgrp s_lgdp96l s_deml, nolog cl(cow)
logit s_ethons05 SPATCONC POPCLDIST s_lr3 ldistpop mtgrp s_lgdp96l s_deml if CONC!=1, nolog cl(cow)












